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Abstract 

This paper presents a search for Higgs bosons decaying to four leptons, either electrons or 
muons, via one or two light exotic gauge bosons Zj, H ZZ^ ox H ^ Z^Zj ^ 4f. 
The search was performed using pp collision data corresponding to an integrated luminosity 
of about 20 fb”* at the center-of-mass energy of = 8 TeV recorded with the ATLAS 
detector at the Large Hadron Collider. The observed data are well described by the Standard 
Model prediction. Upper bounds on the branching ratio oi H ^ ZZd M and on the 
kinetic mixing parameter between the Z^ and the Standard Model hypercharge gauge boson 
are set in the range (1-9) X 10”^ and (4-17) x 10”^ respectively, at 95% confidence level 
assuming the Standard Model branching ratio of 7/ ^ ZZ* ^ 47, for Zd masses between 
15 and 55 GeV. Upper bounds on the effective mass mixing parameter between the Z and 
the Zd are also set using the branching ratio limits in the H ZZd 47 search, and are 
in the range (1.5-8.7) x 10”^ for 15 < mzj < 35 GeV. Upper bounds on the branching ratio 
of 7/ —> ZdZd —> 47 and on Higgs portal coupling parameter, controlling the strength of the 
coupling of the Higgs boson to dark vector bosons, are set in the range (2-3) x 10“^ and (1- 
10) X 10”^ respectively, at 95% confidence level assuming the Standard Model Higgs boson 
production cross sections, for Zd masses between 15 and 60 GeV. 
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1 Introduction 


Hidden sector or dark sector states appear in many extensions to the Standard Model (SM) [1-10], to 
provide a candidate for the dark matter in the universe [11] or to explain astrophysical observations of 
positron excesses [12-14]. A hidden or dark sector can be introduced with an additional U{\)d dark gauge 
symmetry [5-10]. 

In this paper, we present model-independent searches for dark sector states. We then interpret the results 
in benchmark models where the dark gauge symmetry is mediated by a dark vector boson Z^. The dark 
sector could couple to the SM through kinetic mixing with the hypercharge gauge boson [15-17]. In 
this hypercharge portal scenario, the kinetic mixing parameter e controls the coupling strength of the dark 
vector boson and SM particles. If, in addition, the U{\)d symmetry is broken by the introduction of a dark 
Higgs boson, then there could also be a mixing between the SM Higgs boson and the dark sector Higgs 
boson [5-10]. In this scenario, the Higgs portal coupling k controls the strength of the Higgs coupling 
to dark vector bosons. The observed Higgs boson would then be the lighter partner of the new Higgs 
doublet, and could also decay via the dark sector. There is an additional Higgs portal scenario where 
there could be a mass-mixing between the SM Z boson and Zd [7, 8]. In this scenario, the dark vector 
boson Zd may couple to the SM Z boson with a coupling proportional to the mass mixing parameter 5. 

The presence of the dark sector could be inferred either from deviations from the SM-predicted rates of 
Drell-Yan (DY) events or from Higgs boson decays through exotic intermediate states. Model-independent 
upper bounds, from electroweak constraints, on the kinetic mixing parameter of e < 0.03 are reported in 
Refs. [5, 18, 19] for dark vector boson masses between 1 GeV and 200 GeV. Upper bounds on the kinetic 
mixing parameter based on searches for dilepton resonances, pp ^ Zd ^ below the Z-boson mass 
are found to be in range of 0.005-0.020 for dark vector boson masses between 20 and 80 GeV [20]. The 
discovery of the Higgs boson [21-23] during Run 1 of the Large Hadron Collider (LHC) [24, 25] opens 
a new and rich experimental program that includes the search for exotic decays H ZZd 4f and 
H ZdZd —> This scenario is not entirely excluded by electroweak constraints [5-10, 18, 20]. The 

H ZZd process probes the parameter space of e and niz^, or 6 and mz^, where mz^ is the mass of 
the dark vector boson, and the H ZdZd process covers the parameter space of k and mz^ [5, 6]. DY 
production, pp ^ Zd ^ offers the most promising discovery potential for dark vector bosons in the 
event of no mixing between the dark Higgs boson and the SM Higgs boson. The H ZZd process 
offers a discovery pofential complemenfary fo fhe DY process for mz^ < mz [5, 20]. Bofh of fhese would 
be needed fo undersfand fhe properfies of fhe dark vecfor boson [5] . If fhe dark Higgs boson mixes wifh 
fhe SM Higgs boson, fhe H ZdZd M process would be imporfanf, probing fhe dark secfor fhrough 
fhe Higgs portal coupling [5, 6]. 

This paper presenfs a search for Higgs bosons decaying fo four lepfons via one or fwo Zd bosons using 
pp collision dafa af = 8 TeV collecfed at the CERN LHC with the ATLAS experiment. The search 
uses a dataset corresponding to an integrated luminosity of 20.7 fb“' with an uncertainty of 3.6% for 
H ZZd 4^ based on the luminosity calibration used in Refs. [26, 27], and 20.3 fb“' with an 
uncertainty of 2.8% for H ZdZd 4^ based on a more recent calibration [28]. Same-flavor decays 
of fhe Z and Zd bosons fo elecfron and muon pairs are considered, giving fhe 4e, 2e2p, and Ap final 
sfafes. Final slates including r lepfons are nol considered in fhe H ZZd A£ and H ZdZd Ai 
decays. In fhe absence of a significanl signal, upper bounds are sef on fhe relafive branching ratios 
BRiH ZZd Af)/BR{H A{) and BR{H ZdZd 4^)/BR(// ^ ZZ* ^ Ai) as functions of 
fhe mass of fhe dark vecfor boson mz^- The branching ratio limits are used to set upper bounds on the 
kinetic mixing, mass mixing, and Higgs boson mixing parameters [5, 6]. The search is restricted to the 
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mass range where the from the decay of the Higgs boson is on-shell, i.e. 15GeV < mz^ < mull, 
where = 125 GeV. Dark vector boson masses below 15 GeV are not considered in the present search. 
Although the low-mass region is theoretically well motivated [7, 8 ] , the high px of the boson relative 
to its mass leads to signatures that are better studied in dedicated searches [29]. 

The paper is organized as follows. The ATLAS detector is briefly described in Sec. 2. The signal and 
background modeling is summarized in Sec. 3. The dataset, triggers, and event reconstruction are presen¬ 
ted in Sec. 4. Detailed descriptions of the searches are given in Sec. 5 and Sec. 6 for ^ ZZ^ —> and 
H ZdZd —> 4^ processes respectively. Finally, the concluding remarks are presented in Sec. 7. 

2 Experimental Setup 


The ATLAS detector [30] covers almost the whole solid angle around the collision point with layers of 
tracking detectors, calorimeters and muon chambers. The ATLAS inner detector (ID) has full coverage' in 
the azimuthal angle ^ and covers the pseudorapidity range \r]\ < 2.5. It consists of a silicon pixel detector, 
a silicon microstrip detector, and a straw-tube tracker that also measures transition radiation for particle 
identification, all immersed in a 2 T axial magnetic field produced by a superconducting solenoid. 

High-granularity liquid-argon (LAr) electromagnetic sampling calorimeters, with excellent energy and 
position resolution, cover the pseudorapidity range I 77 I < 3.2. The hadronic calorimetry in the range 
\r]\ < 1.7 is provided by a scintillator-tile calorimeter, consisting of a large barrel and two smaller extended 
barrel cylinders, one on either side of the central barrel. The LAr endcap (1.5 < I 77 I < 3.2) and forward 
sampling calorimeters (3.1 < It/I < 4.9) provide electromagnetic and hadronic energy measurements. 

The muon spectrometer (MS) measures the deflection of muon trajectories with \r]\ < 2.7 in a toroidal 
magnetic field. Over most of the //-range, precision measurement of the track coordinates in the principal 
bending direction of the magnetic field is provided by monitored drift tubes. Cathode strip chambers 
are used in the innermost layer for 2.0 < |?/| < 2.7. The muon spectrometer is also instrumented with 
dedicated trigger chambers, resistive-plate chambers in the barrel and thin-gap chambers in the end-cap, 
covering |?/| < 2.4. 

The data are collected using an online three-level trigger system [31] that selects events of interest and 
reduces the event rate from several MHz to about 400 Hz for recording and offline processing. 


3 Monte Carlo Simulation 

Samples of Higgs boson production in the gluon fusion (ggF) mode, with H ZZ^ 4f and H 
ZdZd —> 4^, are generated for mu - 125 GeV and 15 < mz^ < 60 GeV (in 5 GeV steps) in Mad- 
Graph5 [32] with CTEQ6L1 [33] parton distribution functions (PDF) using the Hidden Abelian Higgs 
Model (HAHM) as a benchmark signal model [5, 9, 10]. PythiaS [34, 35] and Photos [36-38] are used 
to take into account parton showering, hadronization, and initial- and final-state radiation. 


* ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the center of the detector 
and the z-axis along the beam pipe. The x-axis points from the IP to the center of the LHC ring, and the ^-axis points upward. 
The azimuthal angle (j) is measured around the beam axis, and the polar angle d is measured with respect to the z-axis. ATLAS 
defines transverse energy Et = E sin6, transverse momentum pj = p sin0, and pseudorapidity 7/ = -ln(tan[6/2)]. 
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The background processes considered in the H ZZ^ A£ and H Z^Z^ —> Ai searches follow those 

used in the H ZZ* A£ measurements [39], and consist of: 

• Higgs boson production via the SM ggF, VBF (vector boson fusion), WH, ZH, and ttH processes with 
H ZZ* A{ final states. In the H Z^Z^ At search, these background processes are normal¬ 
ized with the theoretical cross sections, where the Higgs boson production cross sections and decay 
branching ratios, as well as their uncertainties, are taken from Refs. [40, 41]. In the H ZZ^ —> At 
search, the normalization of H ^ At is determined from data. The cross section for the ggF pro¬ 
cess has been calculated to next-to-leading-order (NLO) [42-44] and next-to-next-to-leading-order 
(NNLO) [45-47] in QCD. In addition, QCD soft-gluon resummations calculated in the next-to-next- 
to-leading-logarithmic (NNLL) approximation are applied for the ggF process [48]. NLO electroweak 
(EW) radiative corrections are also applied [49, 50]. These results are compiled in Refs. [51-53] as¬ 
suming factorization between QCD and EW corrections. For the VBF process, full QCD and EW 
corrections up to NEO [54-56] and approximate NNEO QCD [57] corrections are used to calculate 
the cross section. The cross sections for the associated WH and ZH production processes are calcu¬ 
lated at NEO [58] and at NNEO [59] in QCD, and NEO EW radiative corrections are applied [60]. The 
cross section for associated Higgs boson production with a ttpair is calculated at NEO in QCD [61-64]. 
The SM ggE and VBE processes producing H ZZ* At backgrounds are modeled with Powheg, 
Pythia8 and CTIO PDEs [33]. The SM WH, ZH, and ttH processes producing H ZZ* At 
backgrounds are modeled with Pythia8 with CTEQ6E1 PDEs. 

• SM ZZ* production. The rate of this background is estimated using simulation normalized to the SM 
cross section at NLO. The ZZ* —> At background is modeled using simulated samples generated with 
Powheg [65] and PythiaS [35] for qq ZZ*, and gg2ZZ [66] and Jimmy [67] for gg ZZ*, and CTIO 
PDEs for both. 

• Z-i-jets and ti. The rates of these background processes are estimated using data-driven methods. How¬ 

ever Monte Carlo (MC) simulation is used to understand the systematic uncertainty on the data-driven 
techniques. The Z-i-jets production is modeled with up to five parfons using Alpgen [68] and is divided 
info fwo sources: Z-i-lighf-jefs, which includes Zee in fhe massless c-quark approximation and Zhh 
wifh hh from parfon showers; and Zhh using mafrix-elemenf calculations fhaf fake info accounf fhe 
fi-quark mass. The mafehing scheme of mafrix elemenfs and parfon shower evolufion (see Ref. [69] 
and fhe references fherein) is used fo remove any double counfing of idenfical jefs produced via fhe 
mafrix-elemenf calculafion and fhe parfon shower, buf fhis scheme is nof implemenfed for f7-jefs. There¬ 
fore, hh pairs wifh separation A/? = + i^rif > 0.4 befween fhe fi-quarks are faken from fhe 

mafrix-elemenf calculafion, whereas for AR < 0.4 fhe parfon-shower hh pairs are used. Eor compar¬ 
ison befween dafa and simulation, fhe NNEO QCD EEWZ [70, 71] and NEO QCD MCEM [72, 73] 
cross-seefion calculafions are used to normalize the simulations for inclusive Z boson and Zhh produc¬ 
tion, respectively. The tf background is simulated with MC@NEO-4.06 [74] with parton showers and 
underlying-event modeling as implemented in Herwig 6.5.20 [75] and Jimmy. The AUET2C [76] tune 
for the underlying events is used for tt with CTIO PDEs. 

• SM WZ and WW production. The rates of these backgrounds are normalized to theoretical calculations 
at NEO in perturbative QCD [77]. The simulated event samples are produced with Sherpa [78] and 
CTIO PDEs. 

• Backgrounds containing 7/i/r and T, namely ZJjxfj and ZT. These backgrounds are normalized using 
the ATEAS measurements described in Ref. [79]. These processes are modeled with Pythia8 [35] and 
CTEQ6E1 PDEs. 

Differing pileup condifions (mulfiple profon-profon inferaefions in fhe same or neighboring bunch-crossings) 
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as a function of the instantaneous luminosity are taken into account by overlaying simulated minimum- 
bias events generated with PythiaS onto the hard-scattering process and reweighting them according to 
the distribution of the mean number of interactions observed in data. The MC generated samples are pro¬ 
cessed either with a full ATLAS detector simulation [80] based on the GEANT4 program [81] or a fast 
simulation based on the parameterization of the response to the electromagnetic and hadronic showers 
in the ATLAS calorimeters [82] and a detailed simulation of other parts of the detector and the trigger 
system. The results based on the fast simulation are validated against fully simulated samples and the 
difference is found to be negligible. The simulated events are reconstructed and analyzed with the same 
procedure as the data, using the same trigger and event selection criteria. 


4 Event Reconstruction 

A combination of single-lepton and dilepton triggers is used to select the data samples. The single¬ 
electron trigger has a transverse energy (Et) threshold of 25 GeV while the single-muon trigger has a 
transverse momentum (px) threshold of 24 GeV. The dielectron trigger has a threshold of Ex - 12 GeV 
for both electrons. In the case of muons, triggers with symmetric thresholds at px - 13 GeV and 
asymmetric thresholds at 18 GeV and 8 GeV are used. Finally, electron-muon triggers are used with 
electron Ej thresholds of 12 GeV or 24 GeV depending on the electron identification requirement, and a 
muon px threshold of 8 GeV. The trigger efficiency for evenfs passing fhe final selecfion is above 97% [39] 
in each of fhe final sfafes considered. 

Dafa evenfs recorded during periods when significanf porfions of fhe relevanf defector subsysfems were 
nof fully functional are rejecfed. These requiremenfs are applied independenfly of fhe lepton final sfafe. 
Evenfs in a fime window around a noise bursf in fhe calorimefer are removed [83]. Furfher, all friggered 
evenfs are required to confain a reconsfrucfed primary verfex formed from af leasf 3 fracks, each wifh 
Px > 0.4 GeV. 

Election candidates consisf of clusfers of energy deposited in fhe electiomagnefic calorimefer and asso- 
ciafed wifh ID fracks [84]. The clusfers mafched fo fracks are required to satisfy a sef of identification 
criferia such fhaf fhe longifudinal and fransverse shower profiles are consisfenf wifh fhose expecfed from 
electiomagnefic showers. The election fransverse momenfum is compufed from fhe clusfer energy and 
fhe frack direcfion af fhe inferacfion poinf. Selecfed elections musf salisfy Ej > 7 GeV and \r]\ < 2.47. 
Each election musf have a longifudinal impacf paramefer (zo) of less fhan 10 mm wifh respecf to fhe re- 
consfrucfed primary verfex, defined as fhe verfex wifh af leasf fhree associafed fracks for which fhe ^ p^ 
of fhe associafed fracks is fhe highesf. Muon candidafes are formed by mafching reconsfrucfed ID fracks 
wifh eifher complefe or partial fracks reconsfrucfed in fhe muon specfromefer [85]. If a complefe frack 
is presenf, fhe fwo independenf momenfum measuremenfs are combined; ofherwise fhe momenfum is 
measured using fhe ID. The muon reconstiucfion and identification coverage is exfended by using fracks 
reconsfrucfed in fhe forward region (2.5 < \r]\ < 2.7) of fhe MS, which is oufside fhe ID coverage. In 
fhe center of fhe barrel region (Ipl < 0.1), where fhere is no coverage from muon chambers, ID fracks 
wifh Px > 15 GeV are idenfified as muons if fheir calorimetiic energy deposifs are consisfenf wifh a 
minimum ionizing parficle. Only one muon per even! is allowed fo be reconstiucfed in fhe MS only or 
idenfified wifh fhe calorimefer. Selecfed muons musf safisfy px > 6 GeV and Ipl < 2.7. The requiremenf 
on fhe longifudinal impacf paramefer is fhe same as for elections excepf for fhe muons reconsfrucfed in 
fhe forward region wifhouf an ID frack. To rejecf cosmic-ray muons, fhe impacf paramefer in fhe bending 
plane (do) is required fo be wifhin 1 mm of fhe primary vertex. 
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In order to avoid double-counting of leptons, an overlap removal procedure is applied. If two reconstruc¬ 
ted electron candidates share the same ID track or are too close to each other in rj and cp (AR < 0.1), the 
one with the highest transverse energy deposit in the calorimeter is kept. An electron within AR = 0.2 of a 
muon candidate is removed, and a calorimeter-based reconstructed muon within AR = 0.2 of an electron 
is removed. 

Once the leptons have been selected with the aforementioned basic identification and kinematic require¬ 
ments, events with at least four selected leptons are kept. All possible combinations of four leptons 
(quadruplets) containing two same-flavor, opposite-charge sign (SFOS) leptons, are made. The selected 
leptons are ordered by decreasing transverse momentum and the three highest-px leptons should have 
respectively px > 20 GeV, px > 15 GeV and px > 10 GeV. It is then required that one (two) leptons 
match the single-lepton (dilepton) trigger objects. The leptons within each quadruplet are then ordered in 
SFOS pairs, and denoted 1 to 4, indices 1 and 2 being for the first pair, 3 and 4 for the second pair. 


5 H ^ ZZd 4e 

5.1 Search strategy 

The H ZZd 4^ search is conducted with the same sample of selected At events as used in Refs. [26, 
27] with the four-lepton invariant mass requirement of 115 < mi,[ <130 GeV. This collection of events 
is referred to as the At sample. The invariant mass of the opposite-sign, same-flavor pair closest to the 
Z-boson pole mass of 91.2 GeV [ 86 ] is denoted m\ 2 - The invariant mass of the remaining dilepton pair 
is defined as The H ^ At yield, denoted n{H At), is determined by subtracting the relevant 
backgrounds from the At sample as shown in Eq. (1): 


n{H At) =n{At) - n{ZZ*) - n{tt) - n{Z + jets). (1) 

The other backgrounds from WW, ITZ, ZJjifj and ZT are negligible and not considered. 

The search is performed by inspecting the m 34 mass spectrum and testing for a local excess consistent with 
the decay of a narrow Z^ resonance. This is accomplished through a template fit of the m 34 distribution, 
using histogram-based templates of the H ZZ^ —> At signal and backgrounds. The signal template is 
obtained from simulation and is described in Sec. 5.2. The ^234 distributions and the expected normaliz¬ 
ations of the tt and Z-i-jets backgrounds, along with the ^234 distributions of the H —> ZZ* —> At back¬ 
ground, as shown in Fig. 1, are determined as described in Sec. 5.4. The pre-fit signal and H ZZ* At 
background event yields are set equal to the H ^ At observed yield given by Eq. (1). The expected yields 
for the At sample are shown in Table 1 . 

In the absence of any significant local excess, the search can be used to constrain a relative branching 
ratio Rb, defined as: 


Rb = 


BRiH ZZd At) 
BRiH At) 

BRiH ZZd 


At) 


BRiH ZZd 40 -t BRiH ZZ* At) ’ 


( 2 ) 
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Channel 

ZZ* 

tt + Z±jets 

Sum 

Observed 


4n 

3.1 ±0.02 ±0.4 

0.6 ± 0.04 ± 0.2 

3.7 ± 0.04 ± 0.6 

12 

8.3 ± 0.04 ± 0.6 

4e 

1.3 ±0.02 ±0.5 

0.8 ± 0.07 ± 0.4 

2.1 ±0.07 ±0.9 

9 

6.9 ± 0.07 ± 0.9 

2 /i 2 e 

1.4 ±0.01 ±0.3 

1.2 ±0.10 ±0.4 

2.6 ± 0.10 ± 0.6 

7 

4.4 ±0.10 ±0.6 

2 e 2 // 

2.1 ±0.02 ±0.3 

0.6 ± 0.04 ± 0.2 

2.7 ±0.10 ±0.5 

8 

5.3 ± 0.04 ± 0.5 

all 

7.8 ± 0.04 ± 1.2 

3.2 ±0.1 ± 1.0 

11.1 ± 0.1 ± 1.8 

36 

24.9 ±0.1 ± 1.8 


Table 1; The estimated pre-fit background yields of (MC) ZZ*, (data-driven) tt + Z-i-jets, their sum, the observed 
event yield and the estimated pre-fit // —> contribution in the 4£ sample. The H ^ 41 estimate in the last column 
is obtained as the difference between the observed event yield and the sum of the ZZ* and tt + Z-rjets backgrounds. 
The pre-fit H —> ZZ* ^ 4( background and H —» ZZj 4( signal events are normalized to the // —> 4^ observed 
events. The uncertainties are statistical and systematic respectively. The systematic uncertainties are discussed in 
Section 5.5. Uncertainties on the H 4i rates do not include the statistical uncertainty from the observed number. 
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Figure 1; The distribution of the mass of the second lepton pair, 1111 , 4 , of the ^fs= 8 TeV data (filled circles with 
error bars) and the expected (pre-fit) backgrounds. The H —> ZZ* —> 4{ expected (pre-fit) normalization, for a mass 
hypothesis of mH-l25 GeV, is set by subtracting the expected contributions of the ZZ*, Z-rjets and ff backgrounds 
from the total number of observed 4i events in the data. 


where Rb is zero in the Standard Model. A likelihood function (X) is defined as a product of Poisson 
probability densities (!P) in each bin (/) of the ^34 distribution, and is used to obtain a measurement of 
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(3) 


^bins 

1=1 

^bins 

= Y\ X {nf +px nf‘‘) + bi{v)), 

1=1 

where pn is the normalization of the H ZZ* —> Al background (and allowed to float in the fit), p the 
parameter of interest related to the H ZZ^ —> Ai normalization and p x pn the normalization of the 
H ZZd —> At signal. The symbol y represents the systematic uncertainties on the background estimates 
that are treated as nuisance parameters, and A^bins the total number of bins of the ^34 distribution. The 
likelihood to observe the yield in some bin, given the expected yield is then a function of the 
expected yields n{H —> At) of H ^ ZZ^ —> At {pn x p x nf‘') and H —> ZZ* —> At (pn x nf'), and the 
contribution of backgrounds bi{v). 

An upper bound on p is obtained from the binned likelihood fit to the data, and used in Eq. (2) to obtain a 
measurement of Rb, taking into account the detector acceptance (,91) and reconstruction efficiency (e): 


Rb 


px ph X n{H At) 


pxpH X n{H 


P 

p + C 


At) + CxpHX n{H 


At) 


(4) 


where C is the ratio of the products of the acceptances and reconstruction efficiencies mH ^ ZZd At 
and H ZZ* At events: 

C = (5) 

-^zz* X ezz* 

The acceptance is defined as the fraction of generated events that are within a fiducial region. The re¬ 
construction efficiency is defined as the fraction of events within the fiducial region that are reconstructed 
and selected as part of the At signal sample. 


5.2 Signal modeling 

A signal would produce a narrow peak in the m 34 mass spectrum. The width of the m 34 peak for the 
Zd signal is dominated by detector resolution for all Zd masses considered. For the individual decay 
channels and their combination, the resolutions of the distributions are determined from Gaussian 
fits. The 0234 resolutions show a linear trend between mz^, = 15 GeV and mz^, = 55 GeV and vary from 
0.3 GeV to 1.5 GeV respectively, for the combination of all the final states. The resolutions of the 22234 
distributions are smaller than the mass spacing between the generated signal samples (5 GeV), requiring 
an interpolation to probe intermediate values of mz^. Histogram-based templates are used to model the 
Zd signal where no simulation is available; these templates are obtained from morphed signals produced 
with the procedure defined in Ref. [87]. The morphed signal templates are generated with a mass spacing 
of 1 GeV. 

The acceptances and reconstruction efficiencies of the H ZZd At signal and H ZZ* At 
background are used in Eqs. (4) and (5) to obtain the measurement of the relative branching ratio Rb- The 
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Method 

Estimated background 

Afi 

m \2 fit: Z+jets contribution 

2.4 ±0.5 ±0.6 

mi 2 fit: tt contribution 

0.14 ±0.03 ±0.03 

leljA 

m \2 fit: Z+jets contribution 

2.5 ± 0.5 ± 0.6 

mi 2 fit: tt contribution 

0.10 ± 0.02 ± 0.02 

2 ^ 2 e 

tt + e-e"^ relaxed requirements: sum of Z + jets and tt contributions 

5.2 ± 0.4 ± 0.5 

Ae 

tt + e-e"^ relaxed requirements: sum of Z + jets and tt contributions 

3.2 ± 0.5 ± 0.4 


Table 2: Summary of the estimated expected numbers of Z+jets and ff background events for the 20.7 fb ' of data at 
Vi = 8 TeV for the full mass range after kinematic selections, for the H —> ZZj ^ search. The first uncertainty 
is statistical while the second is systematic. The uncertainties are given on the event yields. Approximately 80% of 
the tt and Z+jets backgrounds have m 4 { < 160 GeV. 


acceptances and efficiencies are derived with H ZZ 4 —> A£ and H —> ZZ* —> At MC samples where 
the Higgs boson is produced via ggF. The product of acceptance and reconstruction efficiency for VBF 
differs from ggF by only 1.2% and the contribution of VH and tlH production modes is negligible: the 
products of acceptance and reconstruction efficiency obtained using the ggF production mode are used 
also for VBF, VH and tJH. 


5.3 Event selection 

The Higgs boson candidate is formed by selecting two pairs of SFOS leptons. The value of m \2 is required 
to be between 50 GeV and 106 GeV. The value of is required to be in the range 12 GeV < m 34 < 
115 GeV. The four-lepton invariant mass niAi is required to be in the range 115 < niAt < 130 GeV. After 
applying the selection to the 8 TeV data sample, 36 events are left as shown in Table 1 . The events are 
grouped into four channels based on the flavor of the reconstructed leptons. Events with four electrons 
are in the Ae channel. Events in which the Z boson is reconstructed with electrons, and ^34 is formed 
from muons, are in the lelji channel. Similarly, events in which the Z is reconstructed from muons and 
m-iA is formed from electrons are in the Ijjle channel. Events with four muons are in the 4// channel. 


5.4 Background estimation 

The search is performed using the same background estimation strategy as the H ZZ* At meas¬ 
urements. The expected rates of the ti and Z+jets backgrounds are estimated using data-driven methods 
as described in detail in Refs. [26, 27]. The results of the expected ti and Z+jets background estimations 
from data control regions are summarized in Table 2. In the “m \2 fit method”, the m \2 distribution of ti is 
fitted with a second-order Chebychev polynomial, and the Z+jets component is fitted with a Breit-Wigner 
lineshape convolved with a Crystal Ball resolution function [26]. In the “tt + e-e* relaxed requirements” 
method, a background control region is formed by relaxing the electron selection criteria for electrons 
of the subleading pairs [26]. Since a fit to the data using m 34 background templates is carried out in the 
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Systematic Uncertainties (%) 

Source 

A^i 

4e 

2 /r 2 e 

2 e 2 ju 

Electron Identification 

- 

9.4 

8.7 

2.4 

Electron Energy Scale 

- 

0.4 

- 

0.2 

Muon Identification 

0.8 

- 

0.4 

0.7 

Muon Momentum Scale 

0.2 

- 

0.1 

- 

Euminosity 

3.6 

3.6 

3.6 

3.6 

tt and Z-rjets Normalization 

25.0 

25.0 

25.0 

25.0 

ZZ* (QCD scale) 

5.0 

5.0 

5.0 

5.0 

ZZ* {qqfPT)F and a$) 

4.0 

4.0 

4.0 

4.0 

ZZ* (gg/PDE and aj 

8.0 

8.0 

8.0 

8.0 


Table 3; The relative systematic uncertainties on the event yields in the H —> ZZj 4£ search. 


search, both the distribution in ^34 and normalization of the backgrounds are relevant. For all relevant 
backgrounds {H ZZ* At, ZZ*, tt and Z+jets) the ^34 distribution is obtained from simulation. 


5.5 Systematic uncertainties 

The sources of the systematic uncertainties in the H ZZ^ —> At search are the same as in the 
H ZZ* —> At measurements. Uncertainties on the lepton reconstruction and identification efficiencies, 
as well as on the energy and momentum reconstruction and scale are described in detail in Refs. [26, 27], 
and shown in Table 3. The lepton identification is the dominant contribution to the systematic uncertain¬ 
ties on the ZZ* background. The largest uncertainty in the H ZZ^ At search is the normalization of 
the tt and Z-i-jets backgrounds. Systematic uncertainties related to the determination of selection efficien¬ 
cies of isolation and impact parameters requirements are shown to be negligible in comparison with other 
systematic uncertainties. The uncertainty in luminosity [28] is applied to the ZZ* background normaliz¬ 
ation. The electron energy scale uncertainty is determined from Z ^ ee samples and for energies below 
15 GeV from 7/i/^ ^ ee decays [26, 27]. Final-state QED radiation modeling and background contam¬ 
ination affect the mass scale uncertainty negligibly. The muon momentum scale systematic uncertainty 
is determined from Z ^ jufi samples and from J/if/ —> fijj. as well as T ^ /i/r decays [26, 27] . Theory 
related systematic uncertainties on the Higgs production cross section and branching ratios are discussed 
in Refs. [39-41], but do not apply in this search since the H ^ At normalization is obtained from data. 
Uncertainties on the m 34 shapes arising from theory uncertainties on the PDFs and renormalization and 
factorization scales are found to be negligible. Theory cross-section uncertainties are applied to the ZZ* 
background. Normalization uncertainties are taken into account for the Z-f-jets and tt backgrounds based 
on the data-driven determination of these backgrounds. 


5.6 Results and interpretation 

A profile-likelihood test statistic is used with the CLg modified frequentist formalism [88-91] imple¬ 
mented in the RooStats framework [92] to test whether the data are compatible with the signal-plus- 
background and background-only hypotheses. Separate fits are performed for different mz^ hypotheses 
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from 15 GeV to 55 GeV, with 1 GeV spacing. After scanning the 0234 mass spectrum for an excess con¬ 
sistent with the presence of a ^ ZZ^ ^ signal, no significant deviation from SM expectations is 
observed. 

The asymptotic approximation [90] is used to estimate the expected and observed exclusion limits on p 
for the combination of all the final slates, and the result is shown in Fig. 2. The relative branching ratio 
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Figure 2: The observed (solid line) and median expected (dashed line) 95% confidence level (CL) upper limits on 
the parameter p related to the H —> ZZ 4 ^ At normalization as a function of mzj, for the combination of all four 
channels (4p, 4e, 2p2e, 2e2p). The +lcr and +2cr expected exclusion regions are indicated in green and yellow, 
respectively. 

as a function of mz^ is extracted using Eqs. (2) and (4) where the value of C as a function of 27234 is 
shown in Fig. 3, for the combination of all four final slafes. This is fhen used wifh p lo consfrain fhe value 
of Rb, and fhe resull is shown in Fig. 4 for fhe combination of all four final sfafes. 

The simples! benchmark model adds lo fhe SM Lagrangian [ 6 - 8 , 10] a U{\)d gauge symmelry lhaf in- 
Iroduces fhe dark vector boson Z^. The dark vecfor boson may mix kinefically wifh fhe SM hypercharge 
gauge boson wifh kinetic mixing paramefer e [ 6 , 10]. This enables fhe decay H —> ZZ^ Ihrough fhe Hy¬ 
percharge Porfal. The Z^ is assumed fo be narrow and on-shell. Furlhermore, fhe presenl search assumes 
prompl Zrf decays consislenl wifh currenl bounds onefrom elecfroweak conslrainls [18,19]. The coupling 
of fhe Zd to SM fermions is given in Eq. (47) of Ref. [ 6 ] to be linear in e, so that BR(Zrf —> €£) is independ¬ 
ent of e due to cancellations [ 6 ]. In this model, the H ZZd —> 47" search can be used to constrain the 
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0134 [GeV] 

Figure 3: The ratio C of the products of the acceptances and reconstruction efficiencies in ^ ZZ^ —> At and 
H —> ZZ* ^ At events as a function of mj,n. 


hypercharge kinetic mixing parameter e as follows: the upper limit on Rg shown in Fig. 4 leads to an up¬ 
per limit on BR{H ZZ^ —> At) assuming the SM branching ratio of ^ ZZ* —> At of 1.25 x 10“^ [40, 
41] as shown in Fig. 5. The limit on e can be obtained directly from the BR(// —> ZZ^ —> At) upper 
bounds and by using Table 2 of Ref. [5]. The 95% CL upper bounds on e are shown in Fig. 6 as a function 
of nizj in the case e » /<• where k is the Higgs portal coupling. 

The measurement of the relative branching ratio Rg as shown in Fig. 4 can also be used to constrain the 
mass-mixing parameter of the model described in Refs. [7, 8] where the SM is extended with a dark vector 
boson and another Higgs doublet, and a mass mixing between the dark vector boson and the SM Z boson 
is introduced. This model explores how aU{l)d gauge interaction in the hidden sector may manifest itself 
in the decays of the Higgs boson. The model also assumes that the Z^, being in the hidden sector, does 
not couple directly to any SM particles including the Higgs boson (i.e. the SM particles do not carry dark 
charges). However, particles in the extensions to the SM, such as a second Higgs doublet, may carry dark 
charges allowing for indirect couplings via the Z-Z^ mass mixing. The possibility of mixing between the 
SM Higgs boson with other scalars such the dark sector Higgs boson is not considered for simplicity. 
The Z-Zd mass-mixing scenario also leads to potentially observable H ZZd At decays at the LHC 
even with the total integrated luminosity collected in Run 1. The partial widths of ^ ZZd —> 4^ and 
H ZZd are given in terms of the Z-Zd mass-mixing parameter 6 and mz^ in Eq. (34) of Ref. [8] and 
Eq. (A.4) of Ref. [7] respectively. As a result, using the measurement of the relative branching ratio Rg 
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Figure 4; The 95% CL upper limits on the relative branching ratio, Rb — ^ function of nizj. The 

+ 1(T and ±2cr expected exclusion regions are indicated in green and yellow, respectively. 


described in this paper, one may set upper bounds on the product 6 ^ x BR(Zrf ^ 1 £) as function of mz^ as 
follows. From Eq. (2) and for mz^, < {mh - niz) 


BRiH ZZd 4£) _ Rb 

BRiH ZZ* 4{) ~ (1 - Rb)’ 

_ nH^zzd)^^ 

— -X 

Tsm 


BR(Z* ^ !£) X BR(Zd ^ 1 €) 
BR{H ZZ* 4£) 


( 6 ) 


where Fsm is the total width of the SM Higgs boson and T{H ZZ^) Fsm- From Eqs. (4), (A.3) 
and (A.4) of Ref. [7], T{H ZZd) ~ 6^. It therefore follows from Eq. (6), with the further assumption 
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Figure 5; The 95% CL upper limits on the branching ratio of H —> 4^ as a function of mzj using the 

combined upper limit on Rb and the SM branching ratio of H ^ ZZ* —» 

<s: (m^ - m|) that: 


Rb 

i\-RB) 


fimzj) 


6 ^ X BR(Z^ ^ 1£) X 

BR(Z* ^ 2i) /(mzj 
BR{H ZZ* 40 ^ Fsm ’ 
1 (mj^ - m|)^ 

16 ;r 


(V) 


where v is the vacuum expectation value of the SM Higgs field. The limit is set on the product 6 ^ x 
BR(Z^ —> 2£) since both S and BR(Zd —> 20 are model-dependent: in the case where kinetic mixing 
dominates, BR(Zrf ^ 20 ~ 30% for the model presented in Ref. [6] but it could be smaller when Z- 
Zd mass mixing dominates [8]. In the mzj mass range of 15 GeV to {nin - mz), the upper bounds on 
6^ X BRfZrf —> 20 are in the range ~ (1.5-8.7) x 10~^ as shown in Fig. 7, assuming the same signal 
acceptances shown in Fig. 3. 
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Figure 6: The 95% CL upper limits on the gauge kinetic mixing parameter e as a function of mz^ using the combined 
upper limit on the branching ratio of i/ —> ZZ^ ^ A{ and Table 2 of Ref. [5]. 


6 H ^ ZdZd ^ 

6.1 Search strategy 

H ZdZd —> candidate events are selected as discussed in Sec. 6.2. The Z, Jjip, Y vetoes are applied 
as also discussed in Sec. 6.2. Subsequently, the analysis exploits the small mass difference between the 
two SFOS lepton pairs of the selected quadruplet to perform a counting experiment. After the small mass 
difference requirements between the SFOS lepton pairs, the estimated background contributions, coming 
from H ZZ* At and ZZ* At, are small. These backgrounds are normalized with the theoretical 
calculations of their cross sections. The other backgrounds are found to be negligible. Since there is 
no significant excess, upper bounds on the signal strength, defined as the ratio of the H ZdZd At 
rate normalized to the SM H —> ZZ* —> At expectation are set as a function of the hypothesized mz^. 
In a benchmark model where the SM is extended with a dark vector boson and a dark Higgs boson, the 
measured upper bounds on the signal strength are used to set limits on the branching ratio of H ^ ZdZd 
and on the Higgs boson mixing parameter as a function of mz^ [5, 6]. 
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Figure 7; The 95% CL upper limits on the product of the mass-mixing parameter 6 and the branching ratio of Zj 
decays to two leptons (electrons, or muons), 6^ x BR(Z;/ ^ 2i), as a function of mzj using the combined upper limit 
on the relative branching ratio of // —> ZZ^ and the partial width of H ^ ZZ^ computed in Refs. [7, 8 ]. 


6.2 Event selection 

For the H —> search, unlike in the H ZZ* —> At study [93], there is no distinction between 

a primary pair (on-shell Z) and a secondary pair (off-shell Z), since both Z^ bosons are considered to be 
on-shell. Among all the different quadruplets, only one is selected by minimizing the mass difference 
Am = \m \2 - 0234 ! where mi 2 and m 34 are the invariant masses of the first and second pairs, respect¬ 
ively. The mass difference Am is expected to be minimal for the signal since the two dilepton systems 
should have invariant masses consistent with the same mz^. No requirement is made on Am; it is used 
only to select a unique quadruplet with the smallest Am. Subsequently, isolation and impact parameter 
significance requirements are imposed on the leptons of the selected quadruplet as described in Ref. [39]. 
Figure 8 shows the minimal value of Am for the 2e2fj. final state after the impact parameter significance 
requirements. Similar distributions are found for the 4e and 4fi final states. The dilepton and four-lepton 
invariant mass distributions are shown in Figs. 9 and 10 respectively for mi 2 and m 34 combined. 

For the H Z^Zj —> At search with hypothesized mz^, after the impact parameter significance require¬ 
ments on the selected leptons, four final requirements are applied: 

(1) 115 < m 4 ^ <130 GeV where m 4 ^ is the invariant mass of the four leptons in the quadruplet, consistent 
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Figure 8: Absolute mass difference between the two dilepton pairs, Am = |mi 2 - m 34 | in the 2e2^l channel, for 
mu — 125 GeV. The shaded area shows both the statistical and systematic uncertainties. The bottom plot shows 
the significance of the observed number of events in the data compared to the expected number of events from the 
backgrounds. These distributions are obtained after the impact parameter significance requirements. 


with the mass of the discovered Higgs boson of about 125 GeV [94]. 

(2) Z, 7/(/^, and T vetoes on all SFOS pairs in the selected quadruplet. The Z veto discards the event if 
either of the dilepton invariant masses is consistent with the Z-boson pole mass: |nii 2 - mz\ < 10 GeV 
or |m 34 - mz\ <10 GeV. For the Jjijj and T veto, the dilepton invariant masses are required to be above 
12 GeV. This requirement suppresses backgrounds with Z bosons, Jlijj, and T. 

(3) the loose signal region requirement: nii 2 < m///2 and ^34 < m///2, where mn - 125 GeV. In the 

H —> Z^Zd —> search, the kinematic limit for on-shell Z^ is mz^ < mnl^. 

(4) the tight signal region requirement: \mz^ -^ 12 ! < Sm and \mz^ -ni 34 l < Sm. The optimized values of the 
6 m requirements are 5/3/4.5 GeV for the AelAjillelji final states respectively (the 6 m requirement 
varies with the hypothesized mz^, but the impact of the variation is negligible). This requirement 
suppresses the backgrounds further by restricting the search region to within 6 m of the hypothesized 
mz,. 

These requirements (l)-(4) define the signal region (SR) o^ H ^ Z^Z^ ^ A€ that is dependent on the 
hypothesized mz^, and is essentially background-free, but contains small estimated background contribu¬ 
tions from H —> ZZ* —> At and ZZ* ^ At processes as shown in Sec. 6.5. 
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Figure 9: Dilepton invariant mass, nia = mi 2 or m^ 4 , in the combined 4e + 2e2^ + Afi final state, for nin =125 GeV. 
The shaded area shows both the statistical and systematic uncertainties. The bottom plots show the significance of 
the observed number of events in the data compared to the expected number of events from the backgrounds. These 
distributions are obtained after the impact parameter significance requirement. 


6.3 Background estimation 

For the H ZdZd A£ search, the main background contributions in the signal region come from 
the H ZZ* —> A£ and ZZ* processes. These backgrounds are suppressed by the requirements 

of the tight signal region, as explained in Sec. 6.2. Other backgrounds with smaller contributions come 
from the Z+jets and tt, WW and WZ processes as shown in Fig. 11. The H ZZ* —> At, ZZ* —> At, 
WW and WZ backgrounds are estimated from simulation and normalized with theoretical calculations of 
their cross sections. After applying the tight signal region requirements described in Sec. 6.2, the Z+jets, 
tt and diboson backgrounds are negligible. In the case where the Monte Carlo calculation yields zero 
expected background events in the tight signal region, an upper bound at 68% CL on the expected events 
is estimated using 1.14 events [86], scaled to the data luminosity and normalized to the background cross 
section: 


^background < L X CT X | 1, (8) 

\ Atot / 

where L is the total integrated luminosity, cr the cross section of the background process, and Atot the total 
number of weighted events simulated for the background process. 
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Figure 10: Four-lepton invariant mass, in the combined 4e -H lelji -H Aji final state, for nin - 125 GeV. The shaded 
area shows both the statistical and systematic uncertainties. The bottom plots show the significance of the observed 
number of events in the data compared to the expected number of events from the backgrounds. These distributions 
are obtained after the impact parameter significance requirement. 


To validate the background estimation, a signal depleted control region is defined by reversing the four- 
lepton invariant mass requirement with an < 115 GeV or > 130 GeV requirement. Good 
agreement between expectation and observation is found in this validation control region as shown in 
Fig. 12. 


6.4 Systematic uncertainties 

The systematic uncertainties on the theoretical calculations of the cross sections used in the event selection 
and identification efficiencies are faken info accounf. The effecfs of PDFs, as, and renormalizafion and 
facforizafion scale uncertainties on fhe fofal inclusive cross secfions for fhe Higgs producfion by ggF, VBF, 
VH and ttH are obfained from Refs. [40, 41]. The renormalizafion, facforizafion scales and PDFs and as 
uncertainfies are applied fo fhe ZZ* background esfimafes. The uncerfainfies due fo fhe limited number of 
MC evenfs in fhe tt, Z-i-jels, ZJjijj, ZY and WWjWZ background simulafions are esfimafed as described 
in Sec. 6.3. The luminosify uncertainly [28] is applied fo all signal yields, as well as fo fhe background 
yields fhaf are normalized wifh fheir fheory cross secfions. The defector syslemalic uncerfainfies due 
fo uncertainfies in fhe eleclron and muon idenfificalion efficiencies are esfimafed wifhin fhe accepfance 
of fhe signal region requiremenfs. There are several componenfs fo fhese uncertainfies. For fhe muons. 
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Figure 11: Dilepton invariant mass, mu = m \2 or ^34 after the loose signal region requirements described in Sec. 6.2 
for the 4e, Afi and, 2e2p final states combined, for mn = 125 GeV. The data is represented by the black dots, and 
the backgrounds are represented by the filled histograms. The shaded area shows both the statistical and systematic 
uncertainties. The bottom plots show the significance of the observed data events compared to the expected number 
of events from the backgrounds. The dashed vertical line is the kinematic limit (mi 2 or < 63 GeV) of the loose 
signal region requirements as discussed in Sec. 6.2. 


uncertainties in the reconstruction and identification efficiency, and in the momentum resolution and 
scale, are included. For the electrons, uncertainties in the reconstruction and identification efficiency, 
the isolation and impact parameter significance requirements, the energy scale and energy resolution are 
considered. The systematic uncertainties are summarized in Table 4. 


6.5 Results and interpretation 

Figures 11 and 13 show the distributions of the dilepton invariant mass (for mi 2 and ^34 combined) and 
the absolute mass difference Am = |mi 2 -m 34 | after the loose signal region requirements. Four data events 
pass the loose signal region requirements, one in the Ae channel, two in the 4/i channel and one in the 2e2/i 
channel. Two of these four events pass the tight signal region requirements: the event in the Ae channel 
and one of the events in the 4// channel. The event in the Ae channel has dilepton masses of 21.8 GeV 
and 28.1 GeV as shown in Fig. 11, and is consistent with a mass in the range 23.5 < mz^ < 26.5 GeV. 
For the event in the 4yu channel that passes the tight signal region requirements, the dilepton invariant 
masses are 23.2 GeV and 18.0 GeV as shown in Fig. 11, and they are consistent with a Zd mass in the 
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Figure 12: The minimal absolute mass difference for the 2e2// final state. Events are selected after the impact 
parameter significance requirement and ni/^i i (115,130) GeV. The shaded area shows both the statistical and 
systematic uncertainties. The bottom plots show the significance of the measured number of events in the data 
compared to the estimated number of events from the backgrounds. 


Systematic Uncertainties (%) 

Source 

4// 

4e 

2e2/t 

Electron Identification 

- 

6.7 

3.2 

Electron Energy Scale 

- 

0.8 

0.3 

Muon Identification 

2.6 

- 

1.3 

Muon Momentum Scale 

0.1 

- 

0.1 

Euminosity 

2.8 

2.8 

2.8 

ggF QCD 

7.8 

7.8 

7.8 

ggE PDEs and as 

7.5 

7.5 

7.5 

ZZ* Normalization 

5.0 

5.0 

5.0 


Table 4: The relative systematic uncertainties on the event yields in the H ZjZ^i A£ search. 
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Figure 13: Absolute mass difference. Am — \mi 2 - m 34 | after the loose signal region requirements described in 
Sec. 6.2 for the 4e, and, lel/j. final states combined, for ;«// = 125 GeV. The data is represented by the black 
dots, and the backgrounds are represented by the filled histograms. The shaded area shows both the statistical 
and systematic uncertainties. The bottom plots show the significance of the observed data events compared to the 
expected number of events from the backgrounds. 


range 20.5 < mzj < 21.0 GeV. In the mzj range of 15 to 30 GeV where four data events pass the loose 
signal region requirements, histogram interpolation [87] is used in steps of 0.5 GeV to obtain the signal 
acceptances and efficiencies at the hypothesized mzj. The expected numbers of signal, background and 
data events, after applying the tight signal region requirements, are shown in Table 5. 

For each mz^, in the absence of any significant excess of events consistent with the signal hypothesis, 
the upper limits are computed from a maximum-likelihood fit to the numbers of data and expected signal 
and background events in the tight signal regions, following the CLs modified frequentist formalism [88, 
89] with the profile-likelihood test statistic [90, 91]. The nuisance parameters associated to the systematic 
uncertainties described in Sec. 6.4 are profiled. The parameter of interest in the fit is the signal strength 
jUd defined as the ratio of the H ZdZd 4^ rate relative to the SM H ZZ* A€ rate: 

^ (T X BR(// ^ ^ AQ 

[cr X BRiH ZZ* 4^)]sm ' ^ 

The systematic uncertainties in the electron and muon identification efficiencies, renormalization and fac¬ 
torization scales and PDF are 100% correlated between the signal and backgrounds. Pseudoexperiments 
are used to compute the 95% CL upper bound fid in each of the final states and their combination, and 
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Process 

Ae 

An 

lelfi 

H^ZZ* ^ Ae 
ZZ* Ai 

WIT, WZ 

tt 

Zbb, Z-i-jefs 
ZJjtp and ZY 

(1.5 ±0.3 + 0.2) X 10-^ 
(7.1 ±3.6 ±0.5) X 10-^ 

< 0.7 X 10-2 

< 3.0 X 10-2 

< 0.2 X 10-2 

< 2.3 X 10-3 

(1.0 ±0.3 ±0.3) X 10-2 
(8.4 ± 3.8 ± 0.5) X 10-3 

< 0.7 X 10-2 

< 3.0 X 10-2 

< 0.2 X 10-2 

< 2.3 X 10-3 

(2.9 ± 1.0 ± 2.0) X 10-3 
(9.1 ±3.6 ±0.6) X 10-3 

< 0.7 X 10-2 

< 3.0 X 10-2 

< 0.2 X 10-2 

< 2.3 X 10-3 

Tofal background 
Dafa 

< 5.6 X 10-2 

1 

< 5.9 X 10-2 

0 

< 5.3 X 10-2 

0 

H^ZZ* ^ At 

ZZ* At 

vtvt,vtz 

tt 

Zbb, Z-i-jefs 
ZJ/il/ and ZY 

(1.2 ± 0.3 ± 0.2) X 10-2 
(3.5 ± 2.0 ± 0.2) X 10-3 

< 0.7 X 10-2 

< 3.0 X 10-2 

< 0.2 X 10-2 

< 2.3 X 10-3 

(5.8 ± 2.0 ± 2.0) X 10-3 
(4.1 ±2.7 ±0.2) X 10-3 

< 0.7 X 10-2 

< 3.0 X 10-2 

< 0.2 X 10-2 

< 2.3 X 10-3 

(2.6 ± 1.0 ± 0.2) X 10-3 
(2.0 ± 0.6 ± 0.1) X 10-2 

< 0.7 X 10-2 

< 3.0 X 10-2 

< 0.2 X 10-2 

< 2.3 X 10-3 

Tofal background 
Dafa 

< 5.3 X 10-2 

0 

< 5.1 X 10-2 

1 

< 6.4 X 10-2 

0 


Table 5; The expected and observed numbers of events in the tight signal region of the H —> Z^iZd A€ search 
for each of the three hnal states, for the hypothesized mass niz^ - 25 GeV and 20.5 GeV. Statistical and systematic 
uncertainties are given respectively for the signal and the background expectations. One event in data passes all 
the selections in the 4e channel and is consistent with 23.5 < mzj < 26.5 GeV. One other data event passes all the 
selections in the 4// channel and is consistent with 20.5 < niz^ <21.0 GeV. The H —> ZZ* —> A( numbers are 
summed over the ggF, VBF, ZH, WH and tIH processes. 


for each of the hypothesized mzj. The 95% confidence-level upper hounds on the H —> ^ 

rates are shown in Fig. 14 relative to the SM Higgs hoson process H ZZ* ^ 4^ as a function of the 
hypothesized for the comhination of the three final sfafes 4e, 2e2yu and 4/i. Assuming fhe SM Higgs 
hoson producfion cross section and using BR(// ^ ZZ* 4f’)sM = 1.25 x 10“^ [40, 41], upper bounds 
on fhe branching rafio oiH ^ ZdZd —> 4t can be obfained from Eq. (9), as shown in Fig. 15. 

The simplesf benchmark model is fhe SM plus a dark vecfor boson and a dark Higgs boson as discussed 
in Refs. [6, 10], where fhe branching rafio of Zd —» tt is given as a function of mz^. This can be used 
fo convert fhe measuremenf of fhe upper bound on fhe signal sfrengfh jid info an upper bound on fhe 
branching rafio BR(// —> ZdZd) assuming fhe SM Higgs boson producfion cross secfion. Figure 16 shows 
fhe 95% CF upper limif on fhe branching rafio of // —> ZdZd as a funcfion of mzj using fhe combined 
^d of fhe fhree final sfafes. The weaker bound af higher mzj is due fo fhe facf fhaf fhe branching rafio 
Zd —> drops slighfly af higher mzj [6] as ofher decay channels become accessible. The H ZdZd 
decay can be used fo obfain a mz^-dependenf limif on an Higgs mixing paramefer k' [6]: 


K - KX 


m 


H 


\m 


■ m% 


( 10 ) 


where k is fhe size of fhe Higgs portal coupling and ms is fhe mass of fhe dark Higgs boson. The parfial 
widfh of ^ ZdZd is given in ferms of k [5]. In fhe regime where fhe Higgs mixing paramefer dominafes 
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Figure 14; The 95% confidence level upper bound on the signal strength of H 

4{! in the combined 4e + 2e2/r + A/u final state, for mn =125 GeV. The +lcr and +2cr expected exclusion regions 

are indicated in green and yellow, respectively. 



{k » e), ms > mnll, mz^ < mnll and H Z^Z* —> M is negligible, the only relevant decay is 
H ZdZd- Therefore the partial width Y{H Z^Zj) can be written as: 


T{H 


’Z.dZd) - Tsm 


BR(// ^ ZdZd) 

\ ^ ZdZd) 


( 11 ) 


The Higgs portal coupling parameter k is obtained using Eq. (53) of Ref. [6] or Table 2 of Ref. [5]: 
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Figure 15; The 95% confidence level upper bound on the branching ratio of H ^ ZjZj ^ 4^ as a function of mzj, 
in the combined 4e + 2e2/r + 4fi final state, for mn =125 GeV. The +lcr and +2cr expected exclusion regions are 
indicated in green and yellow, respectively. 



Figure 17 shows the upper bound on the effective Higgs mixing parameter as a function of mz/. for 
< ms < 2mii, this would correspond to an upper bound on the Higgs portal coupling in the range 
K ~ (1-10) X 10-"^. 

An interpretation for H ZdZd is not done in the Z-Zd mass mixing scenario described in Refs. [7, 8 ] 
since in this model the rate ofH^ ZdZd is highly suppressed relative to that ofH ^ ZZd- 


7 Conclusions 

Two searches for an exotic gauge boson Zd that couples to the discovered SM Higgs boson at a mass 
around 125 GeV in four-lepton events are presented, using the ATLAS detector at the LHC. 

The H ZZd 4^ analysis uses the events resulting from Higgs boson decays to four leptons to search 
for an exotic gauge boson Zj, by examining the ^34 mass distribution. The results obtained in this search 
cover the exotic gauge boson mass range of 15 < mz^ < 55 GeV, and are based on proton-proton collisions 
data at = 8 TeV with an integrated luminosity of 20.7 fb“^. Observed and expected exclusion limits 
on the branching ratio oiH ^ ZZd 4^ relative fo H ^ Al are estimated for the combination of all 
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Figure 16; The 95% confidence level upper bound on the branching ratio of // —> in the combined 4e + 2e2fi + 

4ju final state, for nin = 125 GeV. The +lcr and +2cr expected exclusion regions are indicated in green and yellow, 
respectively. 


the final states. For relative branching ratios above 0.4 (observed) and 0.2 (expected), the entire mass 
range of 15 < mzj <55 GeV is excluded at 95% CL. Upper bounds at 95% CL on the branching ratio of 
H —> ZZd 4£ are set in the range (1-9) x 10~^ for 15 < nizj < 55 GeV, assuming the SM branching 
ratio of H ^ ZZ* 4£. 

The H ZdZd —> 4£ search covers the exotic gauge boson mass range from 15 GeV up to the kinematic 
limit of m///2. An integrated luminosity of 20.3 fb“* at 8 TeV is used in this search. One data event is 
observed to pass all the signal region selections in the 4e channel, and has dilepton invariant masses of 
21.8 GeV and 28.1 GeV. This 4e event is consistent with a mass in the range 23.5 < mz^ < 26.5 GeV. 
Another data event is observed to pass all the signal region selections in the 4// channel, and has dilepton 
invariant masses of 23.2 GeV and 18.0 GeVThis 4jU event is consistent with a mass in the range 
20.5 < mz^ <21.0 GeV. In the absence of a significant excess, upper bounds on the signal strength (and 
thus on the cross section times branching ratio) are set for the mass range of 15 < mz^ < 60 GeV using 
the combined 4e, 2e2^, 4/r final states. 

Using a simplified model where the SM is extended with the addition of an exotic gauge boson and a dark 
Higgs boson, upper bounds on the gauge kinetic mixing parameter e (when e » k), are set in the range 
(4-17) X 10“^ at 95% CL, assuming the SM branching ratio ofH^ ZZ* —> 4£, for 15 < mz^ < 55 GeV. 
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Figure 17; The 95% confidence level upper bound on the Higgs mixing parameter k x in^f\in^ - \ as a function of 

nizj, in the combined 4e + 2e2ju + 4ju final state, for nin - 125 GeV. The +lcr and +2cr expected exclusion regions 
are indicated in green and yellow, respectively. 


Assuming the SM Higgs production cross section, upper bounds on the branching ratio of H ^ Z^Z^, 
as well as on the Higgs portal coupling parameter k are set in the range (2-3) x 10“^ and (1-10) x 10“^ 
respectively at 95% CL, for 15 < mz^ < 60 GeV. 

Upper bounds on the effective mass-mixing parameter 6^ x BR(Zrf —> €i), resulting from the U{\)d gauge 
symmetry, are also set using the branching ratio measurements in the H —> ZZ^ —> search, and are in 
the range ( 1.5 - 8.7) x 10 “^ for 15 < niz^ < 35 GeV. 
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